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We report the measurement of A and A yields and inverse slope parameters in d + Au collisions 
at y/SNN = 200 GeV at forward and backward rapidities (y = ± 2.75), using data from the STAR 
forward time projection chambers. The contributions of different processes to baryon transport and 
particle production are probed exploiting the inherent asymmetry of the d + Au system. Compar- 
isons to model calculations show that the baryon transport on the deuteron side is consistent with 
multiple collisions of the deuteron nucleons with gold participants. On the gold side HIJING based 
models without a hadronic rescattering phase do not describe the measured particle yields while 
models that include target remnants or hadronic rescattering do. The multichain model can provide 
a good description of the net baryon density in d + Au collisions at energies currently available at 
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the BNL Relativistic Heavy Ion Collider, and the derived parameters of the model agree with those 
from nuclear collisions at lower energies. 



PACS numbers: 25.75.-q, 25.75.Dw 

I. INTRODUCTION 

The production of strange baryons has been studied ex- 
tensively in heavy-ion collisions at the BNL Relativistic 
Heavy Ion Collider (RHIC) Q. So far, these measure- 
ments have been concentrated around midrapidity. In d 
+ Au collisions, the study of particle production away 
from midrapidity is especially appealing since the inher- 
ent asymmetry of these events allows the probing of dif- 
ferent processes for particle production on the deuteron 
side (the side in the d beam direction, positive rapidity) 
and on the gold side (the side in the Au beam direction, 
negative rapidity) of the reaction. The deuteron side of 
the collision is expected to be dominated by multiple col- 
lisions of the incoming deuteron nucleons with gold par- 
ticipants and corresponding nuclear effects; while on the 
gold side, final state rescattering and soft processes in 
the nuclear breakup might contribute significantly. By 
studying the particle yields and the inverse slope param- 
eters of A and A as a function of centrality, these effects 
are investigated. The centrality dependence of the yields 
can help illuminate the role of nuclear effects in the ob- 
served strangeness enhancement in Au + Au collisions 
II- 

In addition to providing information about strange 
particle production in asymmetric collisions, A and A 
production at forward rapidities in d + Au collisions 
also probes baryon transport and nuclear stopping power. 
The study of nuclear stopping power is a fundamental 
issue in heavy-ion physics since this quantity is re- 
lated to the amount of energy and baryon number that 
get transferred from the beam particles into the reac- 
tion zone. This influences the properties of possible new 
states of matter created in these collisions. For these 
studies, measurements at forward rapidities are crucial, 
since incomplete stopping is expected at collider ener- 
gies [f|. This should lead to significant net baryon den- 
sities near beam rapidity; while at lower energies, the 
bulk of the net baryons are concentrated near midrapid- 
ity. At energies available at the CERN Super Proton 
Synchrotron (SPS), comprehensive studies of the rapid- 
ity distribution of net baryons in asymmetric collision 
systems demonstrate that the rapidity loss in these col- 
lisions depends mainly on the thickness of the nuclear 
target @. Recent theoretical work suggests that nuclear 
stopping may arise from gluon brcmsstrahlung in cold 
nuclear matter 0]. A measurement of the mean rapidity 
loss of baryons in central Au + Au collisions at ^snn = 
200 GeV Q and the comparison to such measurements 
in central Pb + Pb collisions at the SPS @ indicate that 
a saturation of stopping is reached in central nucleus- 
nucleus reactions with respect to the center of mass en- 
ergy at the top RHIC energy Q. This is shown by the 



deviation from a linear scaling of the rapidity loss with 
collision energy observed at RHIC energies. 

This paper presents the measurement of A and A par- 
ticle yields and inverse slope parameters at forward (y 
= 2.75 ± 0.25) and backward rapidity (y — -2.75 ± 
0.25) in d + Au collisions at ^/WWn = 200 GeV. By 
comparing the particleyields to model calculations per- 
formed with AMPT [m EH, EPOS Q3], HIJING [3 
and HIJING/BB [T3, [Tj], information about the mecha- 
nisms for particle production in asymmetric collisions is 
gained. The net A yield (dN/dy(A) - dN/dy(K)) as a 
function of centrality is compared to calculations based 
on the multichain model [l6| that was previously suc- 
cessfully applied to lower energy data. This comparison 
indicates that the baryon rapidity loss in d + Au colli- 
sions at RHIC is consistent with the predictions of the 
multichain model. 



II. EXPERIMENTAL SETUP AND ANALYSIS 
TECHNIQUE 

The data discussed here were taken with the solcnoidal 
tracker (STAR) detector [T| at the RHIC accelerator 
facility. The minimum bias trigger used for the data 
in this analysis required at least one beam momentum 
neutron in the zero degree calorimeter (ZDC) in the Au 
beam direction. This trigger accepts 95±3% of the d 
+ Au hadronic cross section [l8[. The main detectors 
for the present analysis were the two radial-drift for- 
ward time projection chambers (FTPCs) [l9[ that cover 
2.5 < < 4.0 in pseudorapidity on both sides of the 
interaction region. The detector that sits on the side of 
the interaction region the gold particles fly towards, and 
thus intercepts gold fragments, is referred to as FTPC- 
Au (negative rapidity). The detector on the deuteron 
side of the experiment is referred to as FTPC-d (positive 
rapidity). Data from the main TPC [Ij of STAR is used 
to determine the event vertex and to provide a measure of 
the collision centrality, based on the charged particle mul- 
tiplicity at midrapidity. This method of collision central- 
ity determination avoids autocorrelations in the analysis, 
since different detectors are used for the measurements 
discussed here and for the centrality definition. The min- 
imum bias data set is thus divided into three centrality 
bins, as suggested in (2f[. The peripheral bin contains 
the 40% - 100% most peripheral events, the mid-central 
bin contains 20% - 40%, and the central bin contains the 
20% most central events, as denned by the charged par- 
ticle multiplicity. The number of binary collisions, the 
number of d and the number of Au participants for each 
of these bins are determined using Monte Carlo Glauber 
calculations incorporating the Hulthcn wave function of 
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FIG. 1: (Color online) a) Invariant mass distribution of A candidates on the deuteron side with estimated background 
distribution, b) and c) show the background subtracted A and A invariant mass distributions. The widths of the peaks are due 
to the limited momentum resolution of the detectors, and are reproduced by simulations. 



centrality {N par t} {N bin } 

minimum bias 8.3 ± 0.4 7.5 ± 0.4 

central (top 20%) 15.7li; 15.0± ;| 
mid-central (20% - 40%) H-2±i;J 10.6t^ 
peripheral (40% - 100%) 5.1 ± 0.4 4.2 ± 0.4 

centrality {N part ,d) (jjgg 

minimum bias 1.6 6.7 

central (top 20%) 2.0 13.7 

mid-central (20% - 40%) 1.9 9.3 

peripheral (40% - 100%) 1.4 3.7 



TABLE I: Mean number of participants and mean number 
of binary collisions for minimum bias events and the three 
centrality classes, determined by Glauber calculations. Also 
given are the mean number of participants separated for the 
deuteron and the gold nucleus. 



the deuteron [22j . Table[T]summarizes the Glauber calcu- 
lation results. Contrary to the case of collisions of large, 
equal-sized nuclei, in d + Au collisions the mean num- 
ber of binary collisions is smaller than the mean number 
of participants since in most cases each Au nucleon only 
participates in one collision. 

After event selection cuts, which required a recon- 
structed primary event vertex along the beam axis within 
50 cm of the center of the detector system, 10 7 minimum 
bias events were accepted in the data sample. The vertex 
reconstruction efficiency was determined to be 93±1% 
[HI . Since the vertex reconstruction efficiency is a strong 
function of the number of tracks at midrapidity only pe- 
ripheral events are affected by not reconstructed vertices. 
All particle spectra and yields are corrected for trigger 
and vertex finding efficiencies. 

The momentum resolution of the FTPCs was deter- 
mined from simulations to be between about 10% and 



20% for single charged tracks in the momentum and ra- 
pidity range covered by the present analysis. The mo- 
mentum resolution for reconstructed A and A is on the 
order of 20%. The binning in px for the A and A spectra 
presented here is chosen accordingly. 

In the FTPCs, A and A are reconstructed using their 
dominant decay modes A — * pir~ and A — > pTr + , which 
have a branching ratio of 64%. A candidates are identi- 
fied via displaced vertices. The FTPCs measure a max- 
imum of 10 points on a track. Due to the high momen- 
tum of particles in the forward rapidity region and conse- 
quently very similar energy loss dE/dx of different parti- 
cle species in the detector gas, particle identification via 
the specific energy loss is impossible for singly charged 
hadrons. Thus, A candidates are formed from all pairs 
of one positive and one negative track which make up 
a possible decay vertex that is well separated from the 
main event vertex. These A candidates are convention- 
ally referred to as V0 due to their topology and charge. 
In the present analysis, a minimum separation of 20 cm 
is required. This large minimum decay length leads to a 
reduction in the overall reconstruction efficiency which 
is corrected for in the analysis. 

In the case of A the positive track is assumed to be 
a p, while the negative track is assumed to be a ir~ . 
For A the positive track is assumed to be a 7r + , while 
the negative track is assumed to be a p. Since the most 
abundantly produced particle species are pions this lack 
of particle identification introduces a considerable com- 
binatoric background to the measurement. Strict cuts 
on the geometry of the assumed daughter tracks and the 
resulting A candidate efficiently reduce this background. 
The cut with the highest discriminating power was on the 
distance of closest approach (dca) of the decay daughters 
to the primary vertex, which should be relatively small 
for the p candidate and large for the n candidate since 
the heavier decay daughter typically carries most of the 
momentum of the original particle and thus points back 
to the primary vertex while the lighter daughter does 
not. This cut selects track pairs originating from a de- 
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cay vertex well separated from the primary vertex and in 
addition reduces the background from K® — ► ir + ir~ by 
favoring asymmetric decays. Additional cuts with high 
discriminating power were on the dca of the resulting A 
candidate and on the distance of the daughter tracks to 
the reconstructed decay vertex. 

The remaining combinatoric background is determined 
by rotating the positive tracks by 180° with respect to 
the negative tracks in an event and then forming A can- 
didates using the same cuts as on real data. With this 
method the original spatial correlation of tracks is de- 
stroyed, removing real A and A from the sample. The 
combinatoric background, due to random track pairs that 
fulfill all analysis cuts, remains, since this depends on the 
track multiplicity and dca distributions which are pre- 
served in the rotated sample. The subtraction of this 
background results in a 15% to 20% correction. A mech- 
anism leading to the loss of A and A signals is a possible 
misidentification of the charge of one of the two decay 
daughters. The probability increases with the particle 
momentum, so it is much more likely for the proton can- 
didate than for the pion candidate. This effect manifests 
itself in like-sign track pairs that fulfill the geometrical 
requirements for the A reconstruction. The size of this 
effect is determined with these like-sign pairs and is cor- 
rected for in the analysis. It is comparable in size to the 
remaining combinatoric background. 

The major source of background remaining after these 
cuts and corrections is from K® — > 7r + 7r~, where one of 
the two daughter pions is assumed to be a proton. For 
the current analysis a full GEANT detector simulation 
with a HIJING [HI generated distribution, where 
both the transverse momentum and the rapidity spectra 
of the particles were taken from the event generator, was 
used to model this background. On the deuteron side 
it was verified that the K® yield taken from HIJING is 
in good agreement with d + Au data in the kinematic 
region covered by the FTPCs. This was done using the 
same analysis procedure as for A and A, however using 
different cuts to allow for the different kinematics of the 
Kg decay. 

For all VOs that pass the cuts, parameters such as the 
rapidity y, transverse momentum px, and invariant mass 
are calculated. In the present analysis, a rapidity range 
of 2.5 < | y\ < 3.0 was chosen since this range is fully 
within the FTPC acceptance over the transverse momen- 
tum range of 0.5 < pr < 2.0 GeV/c used in the analysis. 

Figure []Ji) shows the invariant mass distribution for A 
candidates on the deuteron side in the pr range from 0.5 
GeV/c to 2.0 GeV/c for d + Au minimum bias events. 
Also shown is the background contribution due to K® 
decays estimated from HIJING events. This background 
is subtracted, resulting in the A and A invariant mass 
distributions shown for the deuteron side in Figure [TJd) 
and for the gold side in Figure [Ip). On the gold side, 
the pt range is restricted to 0.7 GeV/c to 2.0 GeV/c, as 
discussed later in Section HTT1 From gaussian fits to the 
central part of the invariant mass distributions a mass 
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FIG. 2: (Color online) A and A spectra on the deuteron and 
on the gold side in d + Au minimum bias collisions. The 
data points on the gold side are multiplied by 2 for better 
visibility. The statistical errors are smaller than the points 
marking the measurements. The curves show a fit with a 
Boltzmann function in transverse mass to the data points. 

of 1.116 GeV/c 2 was determined for both A and A on 
both sides of the collision, in good agreement with the 
literature value. The width of the mass peak, given by 
the a of the fit, is 24 MeV/c 2 for the deuteron side and 
26 MeV/c 2 for the gold side, driven entirely by the de- 
tector resolution. The reconstructed mass is indepen- 
dent of centrality, but shows a slight pr dependence due 
to the pt dependent detector resolution. The variation 
over the studied transverse momentum range is around 
10 MeV/c 2 , with lower values at low px and higher val- 
ues at high px . The observed invariant mass distributions 
are reproduced by a full GEANT simulation taking into 
account the detector response. 

The raw particle yields are extracted by summing up 
the bin contents of the background-subtracted invariant 
mass distributions from 1.08 GeV/c 2 to 1.24 GeV/c 2 . 
To get from the raw particle yields to corrected yields, 
the acceptance and the efficiency for A and A has to be 
determined. This is done by embedding into real d + 
Au events simulated A decays that were run through a 
GEANT model of the detector and a simulator of the 
FTPC response. The reconstruction efficiency for A and 
A in the range 2.5 < \y\ < 3.0 and 0.5 GeV/c < p T < 
2.0 GeV/c is ~ 6% with a small dependence on pr. This 
number includes the effect of detector acceptance and the 
analysis cuts used. It is dominated by the requirement of 
a well-separated decay vertex. The branching ratio of the 
decay into charged particles is 64%, leading to an overall 
efficiency of ~ 4%. 



III. PARTICLE SPECTRA AND YIELDS 

Transverse momentum spectra for A and A in d + Au 
minimum bias collisions at -^/snn = 200 GeV are shown 
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FIG. 3: (Color online) (a) Comparison of the measured A 
yield with model calculations. 

(b) Comparison of the net A yield with model calculations. 
Statistical errors are shown as vertical error bars, the vertical 
caps show the quadratic sum of statistical and systematic er- 
rors including the overall normalization uncertainty. In both 
panels the target and projectile beam rapidities are indicated 
by arrows. 



in Figure [2] for both sides of the collision. An incorrect 
treatment of defective electronics in FTPC-Au during 
data production led to the inclusion of noisy electronics 
channels in the data analysis, affecting the measurements 
at low px in particular. This is due to an excess of low 
Pt tracks that fulfill the dca cuts for decay daughters, 
which manifests itself in a shifting of the reconstructed 
invariant mass at low px- Thus the region below pt = 0.7 
GeV/c is excluded from the analysis on the Au side. Also 
shown are fits to the data with a Boltzmann distribution 
in transverse mass rriT, 



1 d 2 N 
2ttpt dydpr 



C rriT exp{—rriT /T) 



(1) 



where T is the inverse slope parameter of the spectrum 
and C is the overall normalization. The spectra on both 
collision sides agree well with the assumed exponential 
behavior. 

From the fits the total particle yield in the rapidity 
range 2.5 < \y\ < 3.0 is extrapolated and the inverse 
slope parameters are extracted. The missing low px mea- 
surement in the px spectra on the Au side leads to an ad- 
ditional systematic error both in the yield and the inverse 



slope parameter. The A and A yields as well as the in- 
verse slope parameters are determined for minimum bias 
events and the three individual centrality classes: cen- 
tral (0% - 20%), mid-central (20% - 40%) and peripheral 
(40% - 100%) events. The particle yields are corrected 
for acceptance, efficiency and feed-down from weak de- 
cays of hyperons with higher mass. The feed-down con- 
tribution is taken from HIJING simulations. The fraction 
of detected A and A particles originating from decays of 
higher mass hyperons was determined to be 0.1 ± 0.03. 
This number includes the differences in reconstruction ef- 
ficiencies for primary A (A) and for A (A) from hyperon 
decays due to their displaced production point. As usual 
in heavy-ion collisions no correction is applied for the 
contribution from E° decays. Thus all quoted A yields 
consist of the contributions of primary A and S . Table 
im summarizes the particle yields, while Table Mil shows 
the inverse slope parameters determined from the Boltz- 
mann fits to the spectra as well as the A/ A yield ratio 
determined from the particle yields. Within the assump- 
tion of a Boltzmann distribution the observed slope pa- 
rameters translate into mean transverse momenta (pr) 
between 0.74 GeV/c and 0.82 GcV/c. Within that as- 
sumption the fraction of the total yield that is covered 
by the measurement is ~ 64% on the d side and ~ 43% 
on the Au side. 

The systematic errors quoted for the results include 
several contributions. These contributions are the cut 
parameters, the efficiency determination from embed- 
ding, background normalization and feed down correc- 
tions. The size of each contribution is obtained from the 
size of effects on the results from variations of cuts and 
normalizations and from a comparison of measurements 
in different sub-regions of the detectors. The dominating 
contributions are from uncertainties introduced by the 
cut selection and from the efficiency determination. This 
contribution is up to 12% on the deuteron side and as 
large as 20% on the Au side. Since both the background 
and feed down contributions in the raw signal are rela- 
tively small, the large uncertainties on their normaliza- 
tions do not lead to sizeable systematics on the extracted 
yields. Systematics are evaluated separately for the par- 
ticle yields and the A/A ratio. In the determination of 
the net A yield and of the A/ A ratio a significant frac- 
tion of the systematic effects cancel, leading to smaller 
overall systematics in these quantities compared to the 
A and A yields. Especially in the case of the yields the 
systematics on the Au side are considerably larger than 
on the d side due to the aforementioned problems with 
the treatment of noisy electronics channels. 

As an additional systematic check the fits to the spec- 
tra were also performed with an exponential function in 
rriT of the form 



1 d 2 N 
2ttpt dydpT 



C exp{-m T /T') 



(2) 



where T" is the inverse slope parameter of the exponential 
function and C is the overall normalization. The yields 
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centrality 


A dN/dy 


A dN/dy 


net A dN/dy 


deuteron side (y = 2.75) 


min. bias 
top 20% 
20% - 40% 
40% - 100% 


0.067 ± 0.001 (stat) t° Hg (syst) 
0.106 ± 0.003 (stat) to.ou ( s y st ) 
0.094 ± 0.003 (stat) toTit ( s y st ) 
0.045 ± 0.001 (stat) ±° ™ 7 6 (syst) 


0.035 ± 0.001 (stat) ±S'.ooI ( s yst) 
0.054 ± 0.002 (stat) (syst) 
0.047 ± 0.002 (stat) torn? ( s yst) 
0.025 ± 0.001 (stat) ± 0.004 (syst) 


0.032 ± 0.002 (stat) ± 0.004 (syst) 
0.052 ± 0.004 (stat) ± 0.007 (syst) 
0.047 ± 0.004 (stat) ± 0.006 (syst) 
0.020 ± 0.002 (stat) ± 0.003 (syst) 


gold side (y = -2.75) 


min. bias 
top 20% 
20% - 40% 
40% - 100% 


0.118 ± 0.004 (stat) ±g;Sts ( s y st ) 
0.294 ± 0.017 (stat) ±g;g^ (syst) 
0.163 ± 0.008 (stat) t° °H (syst) 
0.048 ± 0.002 (stat) ± 0.012 (syst) 


0.072 ± 0.002 (stat) ± 0.022 (syst) 
0.176 ± 0.010 (stat) ± 0.054 (syst) 
0.096 ± 0.005 (stat) ± 0.029 (syst) 
0.031 ± 0.002 (stat) ± 0.009 (syst) 


0.046 ± 0.004 (stat) ± 0.006 (syst) 
0.118 ± 0.020 (stat) ± 0.015 (syst) 
0.067 ± 0.009 (stat) ± 0.009 (syst) 
0.017 ± 0.003 (stat) ± 0.003 (syst) 



TABLE II: Corrected yields of A, A and net A on both sides of the collision. In addition to the quoted systematic errors there 
is an overall normalization uncertainty of 10% on the particle yields. 



centrality 



A inverse slope [GeV] 



A inverse slope [GeV] 



A/ A yield ratio 



deuteron side (y = 2.75) 



min. bias 0.209 ± 0.003 (stat) ± 0.009 (syst) 

top 20% 0.221 ± 0.005 (stat) ± 0.010 (syst) 

20% - 40% 0.208 ± 0.005 (stat) ± 0.010 (syst) 

40% - 100% 0.202 ± 0.004 (stat) ± 0.009 (syst) 



0.210 ± 
0.224 ± 
0.213 ± 
0.199 ± 



0.004 (stat) ± 

0.007 (stat) ± 

0.007 (stat) ± 

0.005 (stat) ± 



0.009 (syst) 
0.010 (syst) 
0.010 (syst) 
0.009 (syst) 



0.52 ± 0.02 (stat) ± 0.04 (syst) 
0.51 ± 0.02 (stat) ± 0.05 (syst) 
0.50 ± 0.03 (stat) ± 0.05 (syst) 
0.56 ± 0.03 (stat) ± 0.05 (syst) 



gold side (y = -2.75) 



min. bias 0.219 ± 0.005 (stat) ± 0.013 (syst) 

top 20% 0.217 ± 0.006 (stat) ± 0.013 (syst) 

20% - 40% 0.218 ± 0.007 (stat) ± 0.013 (syst) 

40% - 100% 0.221 ± 0.007 (stat) ± 0.013 (syst) 



0.206 ± 
0.210 ± 
0.204 ± 
0.201 ± 



0.005 (stat) ± 

0.007 (stat) ± 

0.008 (stat) ± 

0.008 (stat) ± 



0.012 (syst) 
0.012 (syst) 
0.012 (syst) 
0.011 (syst) 



0.61 ± 0.03 (stat) ± 0.05 (syst) 
0.60 ± 0.05 (stat) ± 0.05 (syst) 
0.59 ± 0.04 (stat) ± 0.05 (syst) 
0.65 ± 0.05 (stat) ± 0.05 (syst) 



TABLE III: Inverse slope parameters determined from Boltzmann fit in vtit for A and A on both sides of the collision and A/ A 
yield ratios. 



extracted with this exponential function are consistent 
with the yields extracted based on the assumption of a 
Boltzmann distribution. Our systematic errors do not 
include yield variations due to different spectrum func- 
tions. 

In addition to the systematic errors associated with the 
analysis there is a 10% overall normalization uncertainty 
in the particle yields [HI]. This uncertainty is included 
in the systematic errors shown in the model comparison 
figures and in the baryon transport study. 

To investigate the mechanisms that drive particle pro- 
duction at forward and backward rapidity in d + Au col- 
lisions, the measured minimum bias yields are compared 
to a variety of model calculations. Figure [3] shows the 
measured A and net A yield compared to model calcula- 
tions. While the A yield is sensitive to the amount of A- 
A pair production, the net A yield is strongly influenced 
by baryon number transport in the collision. Incoming 
nucleons can be transformed to A via the associated pro- 
duction of strangeness, N + N^N + A + K, leading to 
a strong correlation of net A and net baryon numbers. 



The four models used are based on different principles. 
HIJING [IU treats nucleus- nucleus collisions as a super- 
position of individual nucleon-nucleon collisions with a 
paramctrization of nuclear effects such as shadowing and 
does not include final state effects such as scattering in 
the hadronic phase. The HIJING/BB model [13, El is 
based on HIJING but includes baryon junction inter- 
actions for increased baryon number transport. These 
gluonic structures facilitate baryon number transfer over 
large rapidity intervals (HI . AMPT [ljl [n| is a multi- 
phase model that includes a HIJING-like treatment of 
initial nucleon-nucleon reactions as well as a later rescat- 
tering phase. EPOS is a phenomenological approach 
based on a parton model. It incorporates nuclear effects 
via parton ladders and target and projectile remnants. 
None of the model calculations have been tuned to pro- 
vide agreement with the data presented here, except in 
the case of HIJING/BB. Here the version with a string 
tension constant of Ki = 1.5 GeV/fm is used since this 
showed the best agreement with the A yields on both the 
d and the Au side and thus with the observed A-A pair 
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production. For all model comparisons discussed below 
the systematic errors in the model calculations are not 
considered. 

On the dcuteron side, the HIJING description of 
several consecutive nucleon-nucleon reactions that take 
place as the nucleons of the deuteron pass through the 
gold nucleus is assumed to be appropriate. It is expected 
that all models should give good descriptions of particle 
production on the deuteron side. On the gold side, how- 
ever, effects associated with final state rescattering and 
the breakup of the target nucleus are expected to play a 
major role, so the AMPT and EPOS models are expected 
to show a better performance. 

Comparing the A measurements with the rapidity dis- 
tributions predicted by the models, shown in Figure [3^,), 
it is apparent that indeed the deuteron side is very well 
described, with all models yielding the same result. On 
the gold side, AMPT is below the other three models. It 
agrees best with the data, however, also the two HIJING 
models and EPOS are consistent with the measurement. 
In general, all models used give a fair description of the 
A yield and thus of the A-A pair production in minimum 
bias d + Au collisions. 

Larger differences are seen for the net A yield shown in 
Figure [3)}), which is very sensitive to baryon transport, 
since the A has two quarks in common with a nucleon 
and thus can easily be produced from a nucleon via asso- 
ciated production of strangeness. On the deuteron side, 
HIJING/BB shows the best agreement with the data, 
suggesting multiple nucleon-nucleon collisions with addi- 
tional baryon transport are an appropriate description of 
the deuteron side of the collision. On the gold side signif- 
icant differences between the models are apparent. Nei- 
ther HIJING nor HIJING/BB reproduce the measured 
net A yield at negative rapidity, while AMPT and EPOS 
do. This suggests that target related effects, as imple- 
mented in AMPT and EPOS, have a strong influence on 
strangeness production on the Au side. It appears that 
at least either a final state rescattering phase, as imple- 
mented in AMPT, or the inclusion of target remnants in 
EPOS, is necessary to reproduce the observed net A yield 
on the gold side of the reaction. 

Figure [4] shows the A and net A yields for the three 
separate centrality bins compared to calculations with 
HIJING, HIJING/BB (net A only), AMPT and EPOS. 
The agreement of the models with the A yields on the 
deuteron side observed for the minimum bias dataset 
holds for centrality selected collisions. The evolution of 
the A yield as a function of centrality on the Au side 
exceeds the increase predicted by the HIJING model. 
While the A yield in peripheral events agrees best with 
the AMPT prediction, the central result is consistent 
with all three models. In general, the yield increase on 
the gold side significantly exceeds the yield increase on 
the deuteron side with increasing collision centrality. The 
behavior of the net A yield as a function of centrality is 
consistent with the observations in minimum bias colli- 
sions. While HIJING/BB provides the best match to the 



data on the deuteron side, the gold side is not described 
by the HIJING models. EPOS and AMPT are able to 
describe the centrality evolution of the net A yield on 
the Au side. On the deuteron side, all models indicate 
a transition from large transparency to significant stop- 
ping in the probed centrality range. This behavior will 
be further investigated in Section [Tvl 

The minimum bias A/ A yield ratio together with pre- 
dictions from the four models discussed above is shown 
in Figure El As for the net A yields, all models are close 
to the data on the deuteron side with HIJING/BB show- 
ing the best match. On the gold side AMPT and EPOS, 
which both incorporate nuclear effects, can reproduce the 
measurement. 

An interesting feature of the centrality dependence of 
the A/A ratio, shown in the upper panels of Figure [51 
is that while the measured net A yields change signifi- 
cantly with centrality on both sides of the collision in the 
measured rapidity bins, the ratio stays constant within 
statistical errors. This shows that the increase in baryon 
stopping with collision centrality is not reflected in a de- 
crease of the anti-baryon to baryon ratio, at least not 
in the hyperon sector. While the net A yield, given 
by the difference of A and A yield, is directly linked to 
the amount of baryon number transport in the reaction, 
the anti-particle to particle ratio is influenced by baryon 
transport and A-A pair production. The centrality inde- 
pendence of the ratio suggests that baryon number trans- 
port and pair production increase in a similar way with 
increasing collision centrality and thus with the amount 
of nuclear material traversed by the projectile. This is 
shown in the lower panel of Figure [6] with the net A and 
the A yield as a function of centrality on both collision 
sides. This is in line with previous p + A measurements 
with a proton beam of up to 300 GeV on a fixed target 
which showed very similar A/ A ratios for p + Be and p + 
Pb reactions [24j|. These measurements were performed 
on the projectile (p) side of the collision as a function of 
if, defined as P\\ t A/Pmax, where P||,a is the longitudinal 
component of the A momentum and p max is the maximal 
possible longitudinal momentum of the A in the center of 
mass frame (of a nucleon-nucleon system) . The x p range 
of these measurements was ~ 0.2 to ~ 0.4, compared to 
anif of ~ 0.1 for the d side data presented here. 

From the inverse slope parameters listed in Table IIIII 
it can be seen that the inverse slopes of both A and A 
show a collision side dependent behavior with central- 
ity. Within statistical errors, the inverse slope param- 
eter does not change with the number of collisions on 
the gold side. On the deuteron side, an increase with 
centrality and thus with the number of nucleon-nucleon 
collisions the deuteron constituents participate in is ob- 
served. This effect is attributed to an increase of the 
mean transverse momentum of the particles in each of 
the subsequent collisions the deuteron participants suffer 
on their way through the gold nucleus, and agrees with 
the picture of the deuteron side of the reaction discussed 
above. This observation is in contrast to the behavior of 



9 




FIG. 4: (Color online) Comparison of A and net A yields to model calculations for all three centrality classes. Statistical errors 
are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity 
is indicated by arrows. 
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FIG. 5: (Color online) Minimum bias A/A ratio compared to 
model calculations. On the deuteron side HIJING/BB shows 
the best agreement with the results, while on the Au side only 
AMPT and EPOS give a satisfactory description of the data. 



found for inclusive charged hadron yields in d + Au col- 
lisions at mid-and forward rapidity, where no significant 
sensitivity to various model calculations has been ob- 
served [251 ] . Studies at midrapidity suggest that more 
differential measurements, like the ratio of particle pro- 
duction at backward rapidity to forward rapidity as a 
function of pt, can provide additional information on 
the relative contributions of various physical processes 
to particle production [2(| [2?| . Although such a study 
is beyond the scope of the present work, there is a con- 
sistent picture emerging regarding the model preference 
of the d + Au data at both mid- and forward rapidity. 
Specifically, midrapidity studies do not support models 
based on incoherent initial multiple partonic scattering 
and independent fragmentation, such as HIJING. The 
EPOS model, which provides a good match to the mea- 
surements on A production presented here was also found 
to explain the data at midrapidity in d + Au collisions 
across many observables [111 [28| . 



inclusive charged hadrons in d + Au collisions where an 
increase of the mean pt with centrality was observed on 
the Au side, while no centrality dependence was found on 
the deuteron side [2{|. This suggests a difference in the 
behavior of hyperons and charged hadrons, mainly pions, 
in the dynamical evolution of the nuclear collision. 

From the A and A spectra and yields at forward and 
backward rapidities in d + Au collisions, it is indicated 
that the deuteron side of the collision is dominated by 
multiple consecutive nucleon-nucleon collisions of partic- 
ipants of the incoming deuteron with gold nucleons. On 
the gold side, the HIJING models can not reproduce the 
observed net A production, while models including nu- 
clear effects can. This situation is different from that 



IV. BARYON TRANSPORT AND NUCLEAR 
STOPPING POWER 

The discussions in the previous section can be extended 
to a study of baryon transport in d + Au collisions based 
on comparisons to the Multi-Chain Model (MCM) [Tel ] . 
This model predicts the baryon number transport in nu- 
clear collisions based on simple assumptions. To do this, 
the net A yields presented here have to be related to 
the total number of net baryons in the corresponding 
rapidity ranges. This is done using model calculations 
performed with HIJING/BB On the Au side of the 
collision there are clearly some issues with the descrip- 
tion of the net A yield by the HIJING/BB model, as 
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FIG. 6: A/A ratio and net A and A yields as a function of collision centrality on both the deuteron (left) and the gold side 
(right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold 
side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with 
centrality, shown by the net A yield, is matched by the increase of A- A pair production, thus keeping the A/ A ratio constant 
over a wide centrality range. 



discussed above. Thus only the deuteron side in the ra- 
pidity range from to 4.0 is used to extract the following 
correspondence: net baryons = (10 ± 1) x net A. For 
the purpose of this discussion, this is assumed to be valid 
also on the Au side of the reaction. Since the MCM pre- 
dicts net baryon yields, the model output is scaled by 
this parameter before comparing to the data presented 
here. The model curves are obtained by coupling the 
MCM as described in [l|| with probability distributions 
for the number of binary N + N collisions obtained from 
Glauber calculations. 

In the literature nuclear stopping power is usually de- 
scribed by the mean rapidity loss of incoming baryons in 
the nuclear collisions jj], 

Sy = Vbearn ~ (y) , (3) 

where yt> eam is the beam rapidity and (y) is the mean 
rapidity of the projectile baryons after the collision. In 
the MCM, the mean rapidity loss is related to the single 
phcnomcnological parameter a by 

Sy = (n co u-l)/a+l, (4) 

where n co u is the number of collisions the incoming nu- 
cleon suffers. The distribution of the number of collisions 
and the probabilities for one or two nucleons of the pro- 
jectile participating in the reaction arc determined with 
Glauber calculations using the multiplicity based central- 
ity definitions used in the data analysis. The parameter 
a was originally extracted from an analysis of results on 
p + A — > p+X at 100 GeV fixed-target beam energy [3f3 |. 
with a result of a = 3 ± 1. 

Figure [7] shows the measured net A yields on both sides 
of the collision for all three centrality bins together with 
predictions based on the MCM using a = 2.9. Uncertain- 




FIG. 7: (Color online) Net A dN/dy for central, mid-central 
and peripheral events on both the deuteron and the Au side 
of the collision. The data are compared to calculations of 
the distribution of net baryons obtained with the Multichain 
model [l(J with a — 2.9, scaled by 0.1 to account for the 
conversion from net baryons to net A. An overall scale un- 
certainty of 10% on the model curves from this conversion is 
not shown. See text for details. 



ties of the overall scale of the model curves due to the con- 
version from net baryons to net A are on the order of 10% 
and are not shown here. The value of a = 2.9 adopted for 
the figure is the best fit to the results. Good fits are pro- 
vided in the range of a = 2.9±0.5, ignoring uncertainties 
stemming from the conversion from net baryons to net 
A. The data show good agreement with the MCM inde- 
pendent of collision centrality. The range for the model 
parameter a supported by the data is driven mostly by 
the measurements on the deuteron side. On the gold 
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side the net baryon yield is dominated by baryons trans- 
ported from the target rapidity. The rapidity distribution 
of baryons on the Au side is only weakly dependent on 
a, since most participating target (gold) nucleons only 
suffer one single collision in the reaction and thus only 
baryons transported from the projectile side to the target 
side contribute to an a dependence. The model parame- 
ter extracted from the net A data in d + Au collisions at 
^snn = 200 GeV is consistent with that obtained from 
p + A collisions at 100 GeV fixed target energy. 

The good agreement of the MCM with a common pa- 
rameter for reactions with more than an order of magni- 
tude different center of mass energy suggests that the ra- 
pidity loss of the incoming baryons in p(d) + A collisions 
and thus the nuclear stopping power is largely indepen- 
dent of beam energy and to a good approximation only 
a function of the number of collisions over a wide energy 
range. 

In central Au + Au collisions at y/sJm = 200 GeV a 
saturation of the stopping power with energy has been 
observed that was not seen in previous measurements 
at lower energy Q. From the expectations of MCM 
with a = 2.9 it appears that the rapidity loss in cen- 
tral nucleus-nucleus collisions is lower than that in d + 
Au collisions for a comparable number of collisions per 
incoming baryon. An important difference between the 
collisions of large, equal-sized nuclei and collisions of a 
very small nucleus with a large nucleus is that in the lat- 
ter case the nucleons of the small nucleus collide with 
nucleons from the large nucleus that in almost all cases 
have not participated in the reaction before. This is not 
true in the first case, which is characterized by multiple 
collisions of both projectile and target nucleons. This dif- 
ference can lead to differences in the stopping behavior in 
the reaction and could lead to the different observations 
in the two collision systems. 

V. CONCLUSION 

We have presented measurements of A hyperon pro- 
duction in d + Au collisions at ^/snn = 200 GeV at 



forward (y = 2.75) and backward (y = —2.75) rapid- 
ity. The comparison of minimum bias yields of A and 
net A to a variety of model calculations shows that the 
deuteron side is well described by all models used. On the 
gold side, only AMPT and EPOS are able to explain the 
net A results, suggesting nuclear effects have an influence 
on hyperon production on the gold side of the collision. 
The observed centrality independence of the A/ A ratio 
shows that baryon number transport through associated 
production of strangeness and A-A pair production both 
have a similar dependence on the number of participants 
and the number of binary collisions in the reaction. The 
good agreement of the multichain model with the mea- 
sured net A yields using the same parameter as for lower 
energy data suggests energy independence of the nuclear 
stopping power for p(d) + A collisions over a wide range 
in center of mass energy. 
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